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Future space astrometry missions are planned to measure positions and/or parallaxes of celestial 
objects with an accuracy of the order of the microarcsecond. At such a level of accuracy, it will be 
indispensable to take into account the influence of the mass multipole structure of the giant planets 
on the bending of light rays. Within the parametrized post-Newtonian formalism, we present an 
algorithmic procedure enabling to determine explicitly this influence on a light ray connecting two 
points located at a finite distance. Then we specialize our formulae in the cases where 1) the light 
source is located at space infinity, 2) both the light source and the observer are located at space 
infinity. We examine in detail the cases where the unperturbed ray is in the equatorial plane or in 
a meridian plane. 
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I. INTRODUCTION 



Highly precise astrometry and tests of general relativity in the Solar System will require in the foreseeable future 
the measurements of apparent positions of light sources with an accuracy of the order of one microarcsecond (/ias) or 
better. The Global Astrometric Interferometer for Astrophysics (GAIA) Mission is planned to obtain accuracies in 
the range 4-20 /ias [3, B a nd the Space Interferometer Mission (SIM) is designed to reach a differential accuracy of 
0.6 fias on bright stars |17|. 

This context has stimulated several studies devoted to the bending effects due to the mass multipole moments of 
the bodies of the Solar System 0, 0, H, [13, El) [13, [I3|- Owing to the complexity of the calculations, only the 
contribution due to the quadrupole moment J2 of the deflecting mass has been explicitly determined. For a light ray 
grazing Jupiter, it has been shown that this effect amounts to 240 /ias [lo| and may be proposed as a new test of 
general relativity which could be realized with GAIA [1, [l^]- However, rough estimates show that the Jovian higher 
multipole moments J4 and Jg could produce deflections of the order of 10 /xas and 0.1 /xas, respectively. For this 
reason, we propose here a method enabling to determine the deflection effects of any mass multipole moment of a 
static axisymmetric gravitational field within the parametrized post-Newtonian formalism. This method is based on 
an algorithmic procedure found in [ll[ and (isj . The problem is treated in the general case where both the light source 
and the observer are located at a finite distance. It is worthy of note that being a development of results derived 
in [13], [l3| and fld\ by using the so-called Synge's world function, our procedure avoids any integration of the null 
geodesic equations. 

The paper is organized as follows. In Sec. II we give the notations used in this article. In Sec. HI we recall how in 
a static axisymmetric space-time the vector tangent to a null geodesic can be derived from the time transfer function 
T giving the travel time of a photon between two points. In Sec. IV we determine the contribution of each mass 
multipole J„ to the time transfer function T. In Sec. V we carry out the calculation of the contribution of each mass 
multipole to the direction of a light ray at its point of reception, both the light source and the observer being assumed 
to be located at a finite distance. In Sec. VI we focus on the case where the source of the light ray is located at space 
infinity. Indeed, this assumption is sufficient for the light emitted by stars or by extragalactic objects. In Sec. VII we 
assume that both the source of light and the observer are situated at space infinity. We obtain explicit formulae for 
the contribution of each mass multipole to the light deflection. We carry out a detailed calculation of the influence of 
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J2, J3 and J4. We recover the results on the effects of J2 previously obtained in [T^, [Tl|, [12, [ill and We give 
some concluding remarks in Sec. VIII. 

II. NOTATIONS AND CONVENTIONS 

The Lorentzian metric of space-time is denoted by g. The signature adopted for t; is (+ — ~ ~)- We suppose 
that space-time is covered by a global quasi-Galilean coordinate system (x'') = (x*^, a;), where — ct, t being a time 
coordinate, and x ~ (x'). Greek indices run from to 3, and latin indices run from 1 to 3. 

A bold letter denotes an ordered triple. In order to distinguish the triples built with contravariant components of a 
vector from the ones built with covariant components, we systematically use the notation a = (a^, a^, a^) — (a*) and 
6 = (61, 62, 63) — except for the gradient operator, denoted by V as usual. Let us emphasize that these notations 
are purely conventional and that we shall write for example Z = TV -I- A in order to mean that U = + \i. Such a 
mixing of "contravariant" and "covariant" quantities cannot be ambiguous, since these equalities must be understood 
as equalities of functions, and not as equalities between true four dimensional vectors. Given a = (a'), 6 = (5') and 
c = (ci), we use a.h to denote a*&* and a.c to denote a*Ci (Einstein's convention on repeated indices is used in both 
cases). We put \a\ = (^ya^a^y/s and |b| = {S^^hh^Y'"^. 

G is the Newtonian gravitational constant and c is the speed of light in a vacuum. 

III. TIME TRANSFER FUNCTION AND DIRECTION OF A LIGHT RAY 

We assume that the metric does not depend on x'^ . Let xa — (c^Ai^^a) a-nd xb — {ctBjXs) be two events of 
space-time supposed to be connected by an unique light ray Tab- By convention xa and xb denote the emission 
point and the reception point, respectively. The travel time Ib — tA oi a, photon connecting xa and xs is a function 
of Xa and a;^, so we put 

tB -tA = 'T{xA,XB)- (1) 

We call T the time transfer function (relative to the chosen coordinate system) . Following a theorem shown in , 
the direction tangent to the light ray Tab at point a;^ is defined by the 4- vector having the covariant components 

dT 

{Io)b^I, {li)B = -c-—{xa,xb)- (2) 

OXg 

We suppose that the gravitational field is generated by an isolated axisymmetric body. We are only interested in 
calculating the contributions of the mass multipoles to the bending of light at the order So using a standard 

post-Newtonian gauge \X8\ we may content ourselves with a metric given by 

2W / 1 \ 

500 = 1- — + (3) 
50. = fiV (4) 



2W\ 
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9., = -[^ + l—)5^,+0[-,]. (5) 

where is a potential defined by 



p being the rest mass density of the body and 7 the well-known post-Newtonian parameter which describes the 
three-space curvature. Putting 

Rab^xb-xa], Nab^—^^—^, (7) 

Rab 

the time transfer function T{xa, xb) involved in Eq. ([2]) is then determined by (see, e.g., Ref. [l5l |) 



TixA,XB)^-RAB + ^{l + l)RAB [ W{x^o)W)dX + O ( ^ 

c c Jq \c 



(8) 
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the integral being calculated along the segment of a line of ends xa and x b described by the parametric equations 

x^^a-^{\)^ RabNabX + xa, 0<Z<1. (9) 

IV. MULTIPOLE STRUCTURE OF THE TIME TRANSFER FUNCTION 



Throughout this work, the center of mass O of the body is taken as the origin of the quasi-Cartesian coordinates 
and the axis of symmetry is chosen as the x'^-axis. We put 



r ^\x\, rA^\XA\, rB=\XB \ ■ 



(10) 



We assume that the smallest sphere centered on O and containing the body has a radius equal to the equatorial 
radius of the body and that the segment joining xa and xb is, outside this sphere. At any point x such that r > 
W{x) is then given by the multipole expansion 



W{x) 



GM 



OO 



n=2 



k.x 



(11) 



where k denotes the unit vector along the x'^-axis, the P„ are the Legendre polynomials, M is the mass of the body 
and the coefficients J„ are the mass multipole moments. As a consequence the integral involved in the right-hand 
side (r.h.s.) of Eq. ([5]) may be written as (see Refs. [llj and p^ ) 



/ W{x(^o)W)dX = GM 1-^-Jn 

•^0 L n=2 



gn 



(9x3)" 



F{x,xa,xb) 



01=0 



where F{x, xa, xb) is the Shapiro kernel function defined by 



Fix. XA, xb) = m ( |"-"^lt|"-"^|+;;^^ ) . 

Rab \\x - xa\ + \x - xb\ -Rab J 

Substituting Eq. (fT^ into Eq. ([5]), and then using Eq. yield the expansion 

^ oc 

T{xa,xb) = - Rab +Tm{xa,xb) + Xj„{xa,xb), 



n=2 



where Tm is the well-known Shapiro time delay 



T f \ / (rA+rB±RAB\ 

Tm{xa,xb) = (7+ l)^-ln ' — 

\rA + rB - Rab J 



and each Tj^ is determined by 



^ , . , , ,,GM 1 „ 9" f\x-XA\ + \x-XB\+RAB 

TjJxA, xb) = -(7 + l) — -Jnr, In U , | + | , _ | 



(dx^)' 



(12) 



(13) 



(14) 



(15) 



(16) 



To carry out an explicit calculation of the r.h.s. of Eq. (1161) . let us apply Faa di Bruno's formula ^] giving the nth 
derivative of a composite function f{x) — h[u{x)], namely [19| 



^ ii\i2\...im^. \l\ dx ) \2\ dx 

m— 1 ii ,...,ijri 



1 duV fl d2»A'^ 



1 d^u 
m! dx™ 



(17) 



where denotes the summation over the sets of non negative integers zi, 12, im satisfying the pair of 



equations 



ii + 2^2 + 3^3 + ... + mi„i = n, 
ii + 12 + ■■■ + im — n — m + 1, 



(18) 
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with 1 < m < n. This calculation involves the Ith derivatives of \x — xa\ and \x — xb\ with respect to a-'^ at a; = 0. 
It may be seen that 



1 9'|a; — xa\ 



l_^(-l/2) f k.XA 



(19) 



by comparing the Taylor expansion of |a; — xa — &x\ about x at point a; = with the expansion 



\x — XA — 5x\ = \x — xa\ 



1=0 



k.Sx 

\x - Xa\ 



C 



(-1/2) / k.{x ~ Xa) 



\x~xa\ 



(20) 



where bx — 6x'^k and c\ ^^^^ denotes the Gegenbauer polynomial of degree I and of parameter —1/2 (see, e.g., Ref. 
[l|). Using Eq. and the similar expression which would be obtained for the Ith derivatives of \x — xb\, Faa di 
Bruno's formula leads to 



GM 



Tj^{xa,xb) - (7+ 1)— 



1 



1 



{rA + rB - Rab) 



{rA + rB + Rab) 



n— m+1 



Qnrn{xA,XB), (21) 



where Qnrn{xA;XB) is defined by 



enUxA,XB)^i^ir-"' ■^T^^,1[[Si{xa,Xb)T 



with 



.(-1/2) 



Si{xa,xb) = -ttQ 



k.XA 



(-1/2) 



k.XB 

TB 



(22) 



(23) 



An explicit calculation of each Tj^ is easy. Consider, e.g., the case where n = 2. The only sets of non negative 
integers solutions to Eqs. are {ii = 2} for to = 1 and {ii = 0,22 = 1} for m = 2. Then Eqs. (PT|) - (P5)) give 



Tj^{xa,xb) = 



7 + 1 GM J2 r2 i?^^ 



2 r-Ars 1 + riA-riB 



1 — (fe.n^)^ 1 — (fc.n^)^ / 1 1 \ [k.{nA + nB 



rA 



rB 



rA rB J l + UA-riB 



(24) 



We thus recover by a straightforward calculation a formula that we have previously derived from the multipole 
expansion of Synge's world function (see Refs. [l^ and (l6j). 



V. DIRECTION OF A LIGHT RAY AT THE RECEPTION POINT 



For the sake of brevity, let us use the notation 1^ = {{li)B} and put 



XA Xb 

riA = — , riB = . 

rA rB 



Taking into account Eqs. (HH), (HS]) and ([HI Eq. © yields 

Lb = -Nab + M^a, xb) 
where MxatXb) is given by the multipole expansion 



(25) 



(26) 



\{xa,xb) = \m{xa,xb) + Aj^ {xa,xb), 



n=2 



with 



XmK^a.xb) = (7+ 1) -2 



c^TB 1 + riA-nB 



riB - IH Nab 

rA V rA 



(27) 



(28) 
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and 



XjJxa,xb) = (7+1)^- J„r;? <j (n-m+1) 

m— 1 

ub + Nab 



ub - Nai 



{rA +rB - Rai 



\7i — 7ri-\'2 



rA+TB + Rab)''-"'+'^ 
1 



Qnm{xA,XB) 



1 



where 



rrUxA,XB) = {-ir-"' 



[rA + TB - RabY-"'+^ {ta + VB + i?AB)"-'"+l 

' [n — ni) ! 



(29) 



Yil\^l(^A,XB)T ^W[Sq{xA-,XB)X 



[Pi-i {k.UB) k- Pi {k.nB)nB] 



9=1, 

q^l 



(30) 



By convention, J^™i. [Sq{xA,XB)y'' — 1 when to = 1. 

q=il 

Using the solutions to Eqs. (HH) found at the end of Sec. IV in the case where n = 2 it may be seen that Eqs. ([29 
and dSni) lead to 



Xj^{xa,xb) = (7 + l)^^J2r-e< - [k.{nA + nB)Y 



Ub - Nai 



Ub + Nab 



l-ik.UAf , l-ik.riB) 



2 1 r 



rA TB 
1 (rA+rB)RAB k.{nA + nB, 



Ub - Nai 



[TA+rB-RABY {rA + rB + RabT 
Ub + Nab 



[ta+tb - RabY [ta + tb + Rai 



[k - {k.nB)nB] 



1 RAB'^{k-nB)k+[l-2,{k.nBY]nB 
2^ 



(31) 



r% r\ {I + ua-UbY ' ta I + ua-Ub 

A tedious but straightforward calculation shows that this formula is equivalent to the r.h.s. of Eq. (100) in Ref. 



VI. SOURCE LOCATED AT INFINITY 



Given a unit vector AT, let V be the straight line parallel to N passing through xb- Suppose that point xa is 
moved away from Xb along T> so that lim^^^oo "f^A = N. This limit corresponds to a source at infinity observed at 
Xb in the direction N within the zeroth-order approximation. Then X{xa,xb) becomes a function of TV and xb- So 
we put 

X{N,xb)^ lim X{xa,xb). (32) 

Let 

cosa = —nB-N, < a < tt (33) 

and denote by rc the quantity 

— tb sin a, (34) 

that is the impact parameter of the light ray at the zeroth-order approximation. Moreover, define the unit vector pb 
as 

PB = (2sin|)~'(nB + 7V). (35) 
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Taking the limit of Eqs. ((28)) and ([29]) . and noting that hm^ 
we get 



^{ta+tb — Rab) = fc tanQ;/2 and that Nab = 



2GM 



cos - pB + ^Jn ( — ] cos"+^ - A„(iV,a;B) 

n=2 V'c/ 



(36) 



where 



A2(iV, xb) = [1 - (fc.ns)^ - 4(fc.pB)2] pb + 2{k.pB) [k - {k.nB)nB] 
+ sin| {2(fe.nB)A;+ [l-'S{k.nBf] ub} 



(37) 



and 



AniN,XB) - 2«-2(fc.p^)"-2 {(„_!)[!_ (fc.nB)2] - 4(fc.ps)'} Ps 

-2"-3(fe.ps)«-3 |(„ - 2)[1 - (fc.ns)"] - 4(fc.ps)2} [fc - ik.nB)nB] 
+2"^2 sin |(fe.pB)""' {2(fc.nB)fc + [l - 3(fc.nB)2] ub} 



$rim(n, «m) 



n [c/ '/'^(fc.ns)]" i (n - m + l)(fc.pB)*^PB - |(fe.pB)'^-i [k - {k.nB)nB] 



i 1 

+ sin|(fe.pB)*^5I*' [^^'^'^('^■"b)]"" n [^g"'^'Hfc-«B)]'' [Pz-i(fc.nB)fc-Pi(fe-«B)nB] 



i=2 



9 = 2, 
q^l 



(38) 



for n > 3, the coefficients $r!m(*i, •■•j*t?x) being defined by 

Z Z^^T^), . 



i3 + --- + (m-2)i,^ ^ 



(39) 



The contributions of the multipole moments to the deflection of fight wifi be completely negligible in missions 
like GAIA or SIM except for light rays passing quite nearly the deflecting body. This implies that one can put 
cos a/2 = 1 and neglect the terms involving sin a/2 in the range where the multipole expansion yielded by Eqs. (j36p . 
P7p and p8p is relevant. To justify this approximation, we may content ourselves with discussing Eq. (j37p since 
the quadrupole moment of the giant planets prevails over the higher multipole moments. Considering the case of 
Jupiter, we have — 7.149 x 10^ km and J2 — 0.014736 Q, which implies 4:GM/c^reJ2 = 240 /las. As a consequence 
4:GM/c^rcJ2{i'e/'rc)^ > 1 /las if and only if (iff) Vc < 6.2re. This last inequality implies that the infiuence of J2 may be 
taken into account iff sin a/2 < 5 x 10~^ since the distance tb between Jupiter and GAIA or SIM is always > 6 x 10® 
km. Then |1 — cosa/2| < 2.5 x 10^ and the contribution of the term containing sin a/2 in the r.h.s. of Eq. ([37|) 
is markedly less than one /las. A similar conclusion can be drawn for n > 3. As a consequence, in GAIA or SIM 
missions, the formulae obtained in this section yield results which do not significantly differ from the approximation 
obtained by assuming that the observer is at infinity, as we shall see below. 



VII. SOURCE AND OBSERVER LOCATED AT INFINITY 



Let P be the foot of the perpendicular drawn to V from O. Since rc — \xp\, we may define the unit vector p as 



It is easily seen that pb is given by 



P=— . 40 



a a 

ps = cos — p + sm — AT. (41) 



It follows from Eq. that the case where point B is at infinity on (P) corresponds to a = 0. So, we have 

linVa-^oo Pb — P and we can put 

A(iV,p,rc)- lim X{N,xb). (42) 

rB^oo 
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The corresponding limit of each term A„(Ar, xb) is obtained by replacing a by 0, ub by —TV and pb by p in Eqs. 
([57)1 and ([55]) . It follows from ([55]) that only the solutions to Eqs. (|18p such that + 2i4 + ... + {m — 2)1^ = have 

to be retained in the sum . A simple calculation leads to the multipole expansion 



9 CM 

A(Ar,p,r,) = (7 + l)^ 

where 



P+J2jn[-) A„(Ar,p,r, 



(43) 



A„(Ar,p, r,) = ^ (-1)" , „ [1 - (fc.iV)^] {2(n - m + l)(fc.p)"-2™+2p 

n — 2m + 2)!(to — ll! L 
m=l ^ / \ / 

-(n - 2m + 2)(fc.p)"-2™+i [fe _ {k.N)N] | (44) 

for any n > 2, pn being the integer defined by 

p„ = — + 1 II n IS even, p„ = — ^ — it n is odd. (45) 

In order to discuss the bending of light rays, it is convenient to introduce the orthonormal triad formed by p, N 
and 

q = px N. (46) 

Noting that k - {k.N)N = {k.p)p + {k.q)q and that 1 - {k.Nf ^ {k.pf + (fc.q)^, we find that Eq. (gH) may be 
written as 

A„(iV,p,r,) = i:(-l)- , 'VT^w""'^^. [(^-P)' + ^^■^)T'" h(fc.p)"-^"+^p 
•■^ — ' n — 2m + 2 m — 1 

^{n - 2m + 2)(fc.p)"-2m+i(fc.q)q] . (47) 
Equation (gT]) yields for n = 2, 3, 4: 

A2(Ar,p, r,) = [(fc.q)2 - (fc.p)2] p + 2{k.p){k.q)q, (48) 
A3(Ar,p, r,) = (fc.p) [Sik.qf - [k.pf] p + (k.q) [i{k.pf - [k.qf] q, (49) 
A4(Ar,p, r,) = [6{k.pf{k.qf - {k.pf - {k.qf] p + 4(fe.p)(fc.q) [{k.pf - {k.qf] q. (50) 

Suppose that the observer is at rest at space infinity. This observer sees the light source in the direction AT + A 
determined by the opposite of the space-like contravariant components of the vector tangent to the light ray. It follows 
from Eq. (pS)) and Nab = that A = A(Af, p,rc). So the effectively observed deflection vector is given by the 
multipole expansion 

00 

A = AM(A^,P,r,) + Aj„(Ar,p,re), (51) 

n=2 

where 

AM(Ar,p,r,) (7+l)±^p, (52) 
A,,„(Ar,p,r,) = (7+l)^J„U A„(iV,p,re). (53) 



It follows from Eqs. (I521)-(l531) and (gT]) that A is orthogonal to AT. 

Substituting for A2 from Eq. ([48]) into Eq. (|53p yields the expression of the deflection vector Ajj obtained in 
previous works (see, e.g., Refs. [3] and [11]). 

If line T> lies in the equatorial plane, then fc.p = and q = ±fc. As a consequence Eq. (|47p gives 



A^,,(Ar,p,r,) = (-l)^-+i(7 + l)^J2. (-) P (54) 



2GM 



2k 
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for the mass multipole of even order 2k and 

A,;,,^,(iV,p,r,)^(-l)'(7 + l)^J2m - k (55) 

for the mass multipole of odd order 2Z + 1. We note that Aj^^ and Aj^^^^ are orthogonal whatever k and L 
If I? is in a meridian plane, then k.q — 0. Pointing out that the relation 

V f-lV"— n = -1 

^ (n- 2m + 2)!(m- 1)! 

m=l ' ; V / 

is valid whatever n > 1 [2^ , it may be seen that Eq. (|T7)l reduces to 

A^„(Ar,p,r,) = -(7+l)^J„p (fc.p)> (56) 

c'^T'c \rcj 

It follows from Eq. ([56| that for a ray propagating in a meridian plane the greatest deflecting effect due to the J„ 
occurs when p — ±fe, i.e. when T) is parallel to the equatorial plane. Equation (|56|) shows also that the mass multipole 
moments have no deflection effect when the direction of emission at infinity is parallel to the axis of symmetry, since 
fe.p = in this case. 

The deflection angle S{N,p,rc) is defined as the angle between N and TV + A(iV,p, Tc). Since is a unit vector 
and AT. A = 0, we have 

\5{N,p,n)\ = \A{N,p,r,)\ + (^) • (57) 
Let us briefly examine the contributions of S{N,p,rc) due to J2, J3 and J4. It follows from Eqs. that 

|A,„(iV,p,r,)| ^ (7 + 1)^J„ f'-f-Y [1 - (fe.Ar)2]"/2 (58) 

c^rc \rcj 

for n = 2,3,4 [2l|. In each of these cases, the highest possible value for \Aj^{N,p,rc)\ is reached when TV is 
orthogonal to A;. As a consequence we have 

|A,;„(TV,p,re)|max = (7 + l)^Jnf-') ■ (59) 

c^rc \rcj 

For Jupiter, .h = 0.014736 (see Sec. VI), J3 = 0.000001 and J4 = -0.000587 ^]. With 7 = 1, the predicted 
deflexions of a grazing ray specifically due to J2, J3 and J4 are then in the range 

< 240 ^as, |?,73| < 0.016 /.tas, < 9.6 ^ias. (60) 

So our formulae yield a rigorous confirmation of previous estimates given in Ref. ,10]. The lowest possible value 
for I Aj^(TV,p, rc)\ is and corresponds to the case where TV = ±fc. 



VIII. CONCLUSION 



This work yields a complete determination at the order 1 /c^ of the bending of light in a static gravitational flcld 
generated by an isolated axisymmetric body. The method developed here could be extended to the contributions of 
the spin multipole moments. The explicit formula ()2ip giving the multipole expansion of the time transfer function 
may also be of interest for the analysis of the frequency transfers between two atomic clocks. 
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